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Theories of photoinduced phase transitions have developed along with the progress in ex- 
perimental studies, especially concerning their nonlinear characters and transition dynamics. 
At an early stage, paths from photoinduced local structural distortions to global ones are ex- 
plained in classical statistical models. Their dynamics are governed by transition probabilities 
and inevitably stochastic, but they were sufficient to describe coarse-grained time evolutions. 
Recently, however, a variety of dynamics including ultrafast ones are observed in different 
electronic states. They are explained in relevant electronic models. In particular, a coherent 
lattice oscillation and coherent motion of a macroscopic domain boundary need appropriate 
interactions among electrons and lattice displacements. Furthermore, some transitions proceed 
almost in one direction, which can be explained by considering relevant electronic processes. 
We describe the history of theories of photoinduced phase transitions and discuss a future 
perspective. 
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1. Introduction 

Nonequilibrium phases can be generated from an equi- 
librium one by external stimulations. Photoirradiation is 
one of them and sometimes allows the emergence of a 
new phase, which cannot be reached by simply changing 
temperature or pressure, because the energy of a photon 
is much higher than thermal energies. Although these 
phases finally relax to the starting equilibrium phase, 
they are locally stable, separated from the initial one by 
a substantial free-energy barrier (except for the transi- 
tion from a Mott insulator to a metal, as will be men- 
tioned later). When the relaxation time is long enough, 
the properties of this nonequilibrium phase can be stud- 
ied by laser spectroscopy techniques, which have recently 
made great progress. Thus, photoinduced phase transi- 
tions have been studied extensively, both experimentally 
and theoretically. 1-3 Here, we review progress in theories 
for such phase transitions. 

Photoirradiation creates electrons and holes, or ex- 
citons if they are bound. In most of the cases, they 
are accompanied by local structural deformation. To de- 
velop it spatially, some cooperativity is needed. One 
would ask first what kind of interaction possesses such 
cooperativity. Section 2 introduces classical models, in 
which photoexcited states proliferate in some conditions. 
Ground and excited states are simply assigned to Ising 
variables, so that the electronic transitions are statisti- 
cally described on the basis of diabatic potentials. Sec- 
tion 3 first introduces adiabatic pictures for electronic 
models, which are useful in discussing necessary condi- 
tions for the spontaneous growth of photoexcited states. 
Transition dynamics do not always proceed adiabatically. 
Next, this section introduces the electronic dynamics de- 
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rived from the time-dependent Schrodinger equation and 
shows their relevance to experimental findings. Section 4 
concludes and presents future issues. 

2. Statistic Theory of Photoinduced Phase 
Transitions 

2.1 Linear chain system 

In order to clarify how the local structural distortions 
lead to global ones, a theory based on stochastic pro- 
cesses was proposed at first. Hanamura and Nagaosa in- 
troduced a simple model: 4 ) 
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where E-pc is the Franck-Condon excitation energy, Qi 
the relevant displacement that dominantly interacts with 
an electron in the photoexcited | ei) state at the Zth 
molecule and is treated as a classical variable, and VSQi 
the Stokes shift in this state. The coupling strength be- 
tween the displacements of the Zth and I'th molecules is 
denoted by Ku> with its diagonal clement Ku set at zero. 
Linear chain systems are first investigated in ref. 4. Sup- 
pose that electrons on m consecutive molecules (e.g., on 
sites 1 to to) are in the excited state, and others are in 
the ground state. Each molecule is located on the respec- 
tive equilibrium position. This state is governed by the 
Hamiltonian Hg . Then, consider that an electron on 
a neighboring molecule (e.g., on site 0) is excited. The 
corresponding molecule is now away from the equilib- 
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Fig. 1. Energy diagrams of Oth molecule, which is closest to clus- 
ter of m excited molecules, in ground and excited states. 4,5 ' 

rium position. This state is governed by the Hamiltonian 
Hc n+1 \ The difference between these Hamiltonians is 
given by 

H (m+1) = H (m) + Eyc _ ^g Qo (3) 

Then, an interaction mode go (on site 0) is introduced, 
which is a linear combination of the displacements in 
the state with m excited molecules. The diabatic poten- 
tials for H g "^ and i?| m+1 ^ are drawn as a function of qo 
(Fig. 1). They are displaced parabolas. The relative sta- 
bility is determined by the minima of the two parabolas 
E s = H g m) (q = q g ) and E c = H^ m+1) (q = q c ), where 
q e and q c give the respective minimum points, and by 
ee H g m \q = q c ) and £ c v = H^ m+1 \q a = q g ). If the 
inequality E c > E^ (E g > E^) is satisfied, the excited 
(ground) state at the Oth molecule is unstable and would 
decay into the ground ("excited") state. Otherwise, the 
smaller one of E g and E c corresponds to the stable state, 
and the larger one to the metastable state. The analysis 
based on these relations reveals how a cluster of pho- 
toexcited molecules is stabilized by the electron-lattice 
interaction, whose dimcnsionlcss strength is the lattice 
relaxation energy S divided by Epc > and by the coupling 
K among the displacements. A similar analysis also ex- 
plains how a molecular excitation far from the cluster is 
attracted to the cluster. 

2.2 Mean-field theory 

Nagaosa and Ogawa have extended this theory to gen- 
eral dimensions and considered the optical and ther- 
mal transitions. 5 - 1 When the relaxation of the interaction 
mode Qi is rapid, the system is always in equilibrium 
with a given electronic state. These modes can be inte- 
grated out to reduce the original model (1) to the Ising 
model, where the spin-up and spin-down states corre- 
spond to the excited and ground states: 

H = -hj^vi - \Y. JwW > ( 4 ) 

l tyl' 

where the field h, the exchange </;;<, and the phonon dis- 
persion LOk are given by 

S 
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with R; being the lattice vector of the Ztli molecule, and 
lo\ = q = 1 — K. The intersite couplings are treated in 
the mean-field approximation. For instance, the thermal 
transition probabilities were assumed to be the prod- 
uct of the attempt frequency and the activation fac- 
tor, P*t g = r- 1 cxp[-(£; c -i; c )/(fcBT)] and P^ c = 
r -1 exp [— (E c — £ , g )/(fc B T)], where E c is the energy of 
the intersection of the two parabolas. The stochastic dy- 
namics of the system under and after the optical pumping 
are investigated on the basis of the rate equation for the 
density of the excited states n e (t), which is a function of 
time t. When stable and metastable states are present, 
it is written as 

±n e (t) = I{t) + - n c (t)} P^ g n e (t) , (8) 

where I(t) is the optical pumping term, being positive 
to create the e state or negative to create the g state. 
The activation energies in the exponents of the thermal 
transition probabilities are quadratic functions of n c (t) 
in the mean-field approximation. The time evolution of 
n c (t) after n e (0)=0 depends on the strength Iq and du- 
ration to of the optical pumping. When io is less than 
a threshold value ith, however long the optical pumping 
is applied, n c (t) does not exceed a certain value. After 
the pumping is switched off, the system always relaxes 
to the initial state n e (0)=0. When Iq is greater than 7 t h, 
on the other hand, n e (t) monotonically increases without 
saturation for < t < t . After the pumping is switched 
off, n e (t) relaxes to the final value nf—0 for t < t CI and 
to Uf—1 for to > t CI . The critical value of the duration 
time t cr , above which the system is switched from rii=Q 
to n/=l, is inversely proportional to Iq for large Iq and 
diverges as i cr ~ (Iq — ith) near the diverging point 
Ith ■ The relaxation times are also discussed. The theory 
is applied to photoinduced structural transformations in 
polydiacetylenes. 

2.3 Other theories 

The theories above are applied to systems with weakly 
interacting components to which thermal fluctuations 
are so dominant that the transition probability is a key 
quantity. Spin-crossover complexes are such examples, 
where spins indirectly interact with each other through 
couplings with lattice displacements. The Hanamura- 
Nagaosa theory is extended to systems with internal de- 
grees of freedom, e.g., spin-crossover complexes that are 
composed of dimers and show two-step transitions. 6 ' In 
general, stochastic dynamics can be treated in classical 
models with master equations 7 ~ 10 - ) and Monte Carlo sim- 
ulations. 7 ' 10 - 12) 
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Fig. 2. (Color online) Schematic structural transformation. The 
standing and falling dominoes correspond to the displacements 
accompanied with stable and metastable localized electronic 
states. 



3. Electronic Theories of Photoinduced Phase 
Transitions 

3.1 Domino effect 

Photoinduced dynamics are sometimes regarded as 
domino effects. This description is inappropriate for sys- 
tems showing coherent dynamics, but gives a hint for the 
mechanisms of photoinduced phase transitions. Koshino 
and Ogawa theoretically show a domino-like structural 
transformation in a one-dimensional model composed of 
localized electrons and classical lattice displacements: 13 ' 
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where Uj and duj/dt denote the jth displacement and 
velocity, respectively, fey the elastic constants, e the en- 
ergy difference between the two electronic states at Uj=0, 
to the overlap integral between them, and 7 the electron- 
lattice coupling constant. After a photon is absorbed by 
an electron at a site, the lattice locally relaxes under fric- 
tion to the minimum on the adiabatic potential. Then, 
after a photon is spontaneously emitted, three evolution 
patterns are shown to exist. When the elastic couplings 
kij are too weak, the local structural distortion remains 
locally. When they are too strong, the initial local distor- 
tion is pulled back by surrounding lattice displacements 
to the position before the absorption. Only when they are 
intermediately strong and short-ranged will the initial lo- 
cal distortion trigger a global structural transformation 
as the domino effect (Fig. 2). 

3.2 Adiabatic relaxation path 

Photoinduced phase transitions are often related to 
multistability, where different electronic phases are sta- 
ble and metastable. Such stability may be manifested by 
the presence of a first-order phase transition induced by 
changing temperature or pressure. The metastable phase 
may be hidden in equilibrium so that it is manifested only 
by photoirradiation, which gives a much higher energy 
to the system than thermal energies. The photoinduced 
phase transition dynamics depend on the electronic state 
and especially on how the initial state is prepared. For 
this purpose, relevant itinerant-electron models are nec- 
essary. They have finite transfer integrals between neigh- 
boring molecular orbitals. They are off-diagonal elements 
of the model Hamiltonian, giving transition amplitudes. 
In contrast, the stochastic dynamics within classical sta- 
tistical models are governed by transition probabilities, 



which are often given by Boltzmann factors satisfying 
detailed balance. 

Organic molecular tetrathiafulvalene-chloranil (TTF- 
CA) crystals are the most intensively studied. The donor 
TTF and acceptor CA molecules are alternately stacked 
along the most conducting axis. At low temperature or 
under high pressure with contraction, these molecules be- 
come ionic due to the long-range Coulomb interaction. 14 ) 
Otherwise, they are neutral due to the large difference be- 
tween the ionization potential of the donor molecule and 
the electron affinity of the acceptor molecule. In the ionic 
phase at ambient pressure, these molecules are dimerized. 
Both ionic-to-neutral and neutral-to-ionic transitions are 
photoinduced, as shown by the optical reflectivity. 15 ' It 
is found that photons with 2.3 eV corresponding to in- 
tramolecular excitations at TTF sites are much more 
efficient for generation of neutral domains in the ionic 
background than those with 0.6 eV corresponding to in- 
trachain charge-transfer excitations. 1 6 ' The latter clearly 
have a threshold in the photoexcitation density, below 
which a macroscopic neutral domain cannot be gener- 
ated. Thus, the photoinduced dynamics are very sensi- 
tive to the initial condition of the electronic state, from 
which the lattice relaxation starts. 

These observations have motivated theoretical stud- 
ies based on itinerant-electron models. Huai et al. have 
clarified the adiabatic relaxation path from the Franck- 
Condon-type photoexcited state above the ionic ground 
state to the formation of a large neutral domain, us- 
ing a one-dimensional extended Peierls-Hubbard model 
with alternating potentials (and the interchain interac- 
tion mentioned later): 17 ' 



H 



h.c 



+ U ni^ni,i + Vj (qi,qi +1 )6ni6ni +1 
1 1 

+ Ey(«-9 + i) 2 + Ej(«-« + i) 4 , (io) 
1 1 

where c, (c; )CT ) is the creation (annihilation) operator 
of an electron with spin a at site I, n; )CT = c la ci t(T , 
ni = + ni,i, Sni = ni — 2 for odd I, Sni — m for 
even /, qi is the dimensionless displacement of the Ith 
molecule along the chain from its equidistant position. 
The distance between the Ith and (I + l)th molecules is 
given by d^+i = d (l + q l+1 - qf), where d is the av- 
erage intermolecular distance. The donor and acceptor 
molecules are located at odd and even sites, respectively. 
The nearest-neighbor repulsion strength is assumed to 
depend nonlinearly on the intermolecular distance as 

Vi{q u qi +1 ) = Vo + fli{qi - qi+i) + P 2 {qi - qi+lf , (11) 

where Vb is for the regular lattice, and f3\ (fa) is the 
linear (quadratic) coefficient. The parameter to denotes 
the transfer integral, A the level difference between the 
neighboring orbitals in the neutral limit, and U the on- 
site repulsion strength. The elastic energy is expanded 
up to the fourth order: the parameters Si and S2 are the 
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Fig. 3. Schematic electronic and lattice structures in (a) neutral 
and (b) ionic phases. 
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Fig. 4. Schematic excited domains in ionic background. 



linear and nonlinear elastic constants. 

In the neutral phase, the orbital of the donor molecule 
is almost doubly occupied, while that of the acceptor 
molecule is almost empty [Fig. 3(a)]. The total charge 
of the donor molecule at site I is +(2 — ni) = —Sni, 
while that of the acceptor molecule is — n; = — Sni. In 
the ionic phase, both orbitals are almost singly occupied 
[Fig. 3(b)]. The values of the eight parameters are so de- 
termined that they reproduce the ah initio estimation 
of the transfer integral, the ionicity in the ionic phase, 
that in the neutral phase, the degree of dimcrization in 
the ionic phase, the energies and the relative strength of 
the charge-transfer absorption peaks in the ionic phase, 
and the charge-transfer absorption energy in the neu- 
tral phase. The lattice relaxation path starting from the 
Franck-Condon state with a single charge-transfer exci- 
ton is assumed to be described by a domain, 



qi = (-!/» HA? 



tanh 



l\- l °< 
1 2 ' 
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. (12) 

where (— l)'<?o denotes the uniform dimerization in the 
ionic state, Aq the amplitude of a local distortion in- 
duced by an excited domain, the width of the domain 
boundary, and Iq the domain size (Fig. 4). An interchain 
interaction between displacements is introduced to con- 
fine the domain in the ionic background: 

E ^(m-i-lYqof , (13) 



inter 



=1,2,3 



where qi denotes the displacement of a central chain sur- 
rounded by ionic environmental chains, and Ki the 2ith 
coefficient with respect to the distortion in the central 
chain. 

The adiabatic energy surface of the ground state and 



that of the first excited state are drawn as a function 
of the amplitude Aq and the size Iq of the domain with 
the width of the domain boundary optimized. 17 ^ The 
energy surface of the ground state has a plateau corre- 
sponding to a neutral domain separated from the ground 
state by a low potential barrier. In the energy surface of 
the first excited state, a local minimum corresponding to 
the charge-transfer exciton appears around the Franck- 
Condon state. Another local minimum appears and cor- 
responds to a neutral domain, which are separated from 
the Franck-Condon state by such a high potential barrier 
that the absorption of a single photon with 0.6 eV can- 
not overcome it. The domain is stable only when its size 
exceeds a critical value because the energy required for 
creation of the two domain boundaries makes a small do- 
main unstable in the excited state. The spin- and charge- 
density distributions of the ground and excited states 
with a neutral domain are very similar: they are different 
only at the domain boundaries. Because a single charge- 
transfer exciton cannot trigger the formation of a macro- 
scopic neutral domain, a large excess energy is needed at 
the very beginning of the relaxation to overcome the high 
barrier, for charge-transfer excitons to proliferate, and fi- 
nally to form a macroscopic neutral domain. The process 
described above is only for the early nucleation stage of 
the photoinduced macroscopic transformation. 

3.3 Transition dynamics in TTF-CA 

Along with the progress in experimental techniques, 
the pump-probe reflection spectroscopy acquires higher 
time resolution. Ultrafast optical switching from the ionic 
to the neutral state is observed. 18 - 1 Neutral strings are 
produced within 2 ps after a resonant excitation of the 
charge transfer band at 0.65 eV polarized along the stack- 
ing axis. Just below the transition temperature, a macro- 
scopic conversion from the ionic to the neutral states is 
achieved within 20 ps. It is accompanied with coherent 
motion of the macroscopic neutral-ionic domain bound- 
ary with a much longer period than that of the optical 
lattice oscillation. 

Theoretically, real-time dynamics need to be cal- 
culated for further understanding. First, results are 
obtained by solving the time-dependent Schrodinger 
equation of an electron-lattice model without electron- 
electron interaction, where two charge-density-wave 
states with opposite polarizations are degenerate. 19 ' The 
lattice part is treated classically. The photoconversion is 
shown to be nonlinear. During the photoconversion pro- 
cess, macroscopic oscillations are observed, which are es- 
sentially different from the linear modes in equilibrium. 
Then, by adding the kinetic energy of the displacements, 



#kin = -jmiii > 



(14) 



with the Ztli molecular mass to,/, to the model (10), 
Miyashita et al. have calculated real-time dynamics 
of ionicity and dimerization within the unrestricted 
Hartree-Fock approximation after the occupancy of or- 
bitals in the ionic ground state is initially changed. 20 - ) 
As the number of orbitals whose occupancy is changed 
increases, i.e., with increasing photoexcitation density, 
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Fig. 5. Conversion from ionic to neutral phases, as a function of 
number of excited electrons in 100-site chain. 20 ' The nonlinearity 
is strengthened with the passing of time after the photoexcita- 
tion, t=1.5/u>, lb/to, 100/ct) and 250/a;, where u) denotes the bare 
phonon frequency. The initial lattice temperature is T=10~ 2 tQ. 
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Fig. 6. Final ionicity as a function of number of absorbed photons 
in 100-site chain, when electric field of strength eE cx tdo=bw and 
of frequency ui cx t=28ui below linear-absorption peak at 32u> is 
applied to ionic phase. 22 ) oj denotes the bare phonon frequency. 
The initial lattice temperature is T=10 — to- 



more neutral domains are created. Above the threshold 
excitation density, the neutral phase is finally achieved. 
After the photoexcitation, ionic domains with wrong po- 
larization generally appear, which reduce the correlation 
length of the staggered lattice displacement. As the de- 
gree of initial lattice disorder increases, more solitons ap- 
pear between the ionic domains with opposite polariza- 
tions, which obstruct the growth of neutral domains and 
slow down the transition. From the Fourier transform 
of the ionicity as a function of time, three characteris- 
tic time scales are observed: rapid charge-transfer exci- 
tations corresponding to the peak energy of a charge- 
transfer exciton in the random-phasc-approximation 
spectra, slow charge-transfer excitations strongly coupled 
with optical lattice oscillations, and much slower excita- 
tions due to the collective motion of neutral-ionic domain 
boundaries. The last corresponds to the experimentally 
observed, coherent motion of the macroscopic neutral- 
ionic domain boundary, 18 - 1 and is shown to be sensitively 
affected by collisions with the solitons. When the frac- 
tion of conversion from the ionic to the neutral phases 
is plotted as a function of the number of excited elec- 
trons, the curve evolves with time (Fig. 5). Immediately 
after the photoexcitation, it is close to a linear function. 
Above the threshold number, the fraction monotonically 
increases with time until the photoinduced phase transi- 
tion is completed. Below the threshold, the fraction de- 
creases and eventually becomes very small. This result is 
consistent with the experimental finding. 18 -* 

More recently, experimental results have been summa- 
rized covering both the ionic-to-neutral and the neutral- 
to-ionic transitions with different excitation energies, 
excitation densities and temperatures. 21 ) The ionic-to- 
neutral transition induced by the resonant excitation of 
the charge-transfer band proceeds with (1) initial for- 
mation of a confined one-dimensional neutral domain, 
(2) multiplication of such domains to semi-macroscopic 
neutral states by 20 ps, and (3) evolution in the direc- 
tion normal to the sample surface. The dynamics of the 
neutral-to-ionic transition that is induced by a similar 
excitation are clearly different from the above. Although 
one-dimensional ionic domains are initially produced by 



lights, they decay within 20 ps even if the excitation den- 
sity is high. The initial conversion fraction is a linear 
function of the excitation density. 

Thus, we are motivated to incorporate a pulse of os- 
cillating classical electric field, 

E(t) = E cxt smu> cxt t , (15) 

with amplitude E cxt and frequency w cx t for < t < 
2nN ext / '^cxt with integer iV cx t, into the Peierls phase 
of the transfer integral. The time-dependent Schrodinger 
equation is solved again for the one-dimensional extended 
Pcicrls-Hubbard model with alternating potentials. Now, 
the frequency, the amplitude and the duration of the 
pulse can be explicitly and independently varied. The 
dynamics of the ionic-to-neutral transition are indeed 
qualitatively different from those of the neutral-to-ionic 
transition. 

3.4 Difference between I-to-N and N-to-I transitions 

When the dimerized ionic phase is photoexcited, the 
threshold behavior is observed again by plotting the fi- 
nal ionicity as a function of the increment of the total 
energy, i.e., as a function of the number of absorbed 
photons (Fig. 6). 22 > Above the threshold photoexcita- 
tion density, the system finally enters a neutral state 
with equidistant molecules. This is the case even after 
the electric field is switched off. Once produced by the 
field, a neutral domain grows spontaneously. Its growth 
cannot be stopped by switching off the field. When com- 
pared with the threshold behavior derived from the rate 
equation, 5 ^ we do not observe a quantity corresponding 
to the threshold intensity Ith- Even if the amplitude of 
the pulse is small (i.e., even if the pulse is weak), a transi- 
tion seems finally achieved by setting the pulse duration 
long enough. The deterministic dynamics are more effi- 
cient than the stochastic ones because of the restricted 
energy dissipation. When compared with the domino ef- 
fect, the motion of the neutral domain is similar in that 
its growth is spontaneous once triggered. However, it is 
quite different from the motion of dominoes under strong 
friction on the adiabatic potential, 13 ) which is far from 
being coherent. 
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Fig. 7. Final ionicity as a function of number of absorbed pho- 
tons in 100-site chain, when electric field of strength eE c ^tdo=2u) 
and of frequency o; cx t=28a; below linear-absorption peak at 30a; 
is applied to neutral phase. 24 ) u> denotes the bare phonon fre- 
quency. The initial lattice temperature is T=10~ 2 t(). 

For a given amount of the increment of the total en- 
ergy, the lattice dynamics relative to the charge dynamics 
in the ionic-to-neutral transition depend on the charac- 
ter of the pulse. When the pulse is strong and short, the 
charge transfer takes place on the same time scale with 
the disappearance of dimerization. When the pulse is 
weak and long, the dimerization-induced ferroelectric po- 
larization is disordered first to restore the inversion sym- 
metry, and then the charge transfer takes place to bring 
the system to a neutral state. The difference between 
the two time scales increases as the pulse is weakened. 
Between these time scales, a paraelectric ionic phase is 
realized. In the case of intramolecular excitations at 1.55 
eV, such a new ionic phase with disordered polarizations 
is suggested to appear by time-resolved X-ray diffrac- 
tion. 23 ) This similarity may be due to the fact that, 
when the pulse is weak, the supplied energy is not di- 
rectly used to transfer charge density along the chain. 
It is furthermore shown theoretically that infrared light 
polarized along the chain can also induce the ionic-to- 
neutral transition if the excitation density exceeds a criti- 
cal value, which is currently beyond the value experimen- 
tally achieved. 

When the neutral phase is photoexcited, the linear 
behavior is observed by plotting the final ionicity as a 
function of the increment of the total energy (Fig. 7). 24 ) 
The neutral-to-ionic transition proceeds in an uncooper- 
ative manner. If the oscillating electric field is turned off 
before the transition is completed, the ionicity remains 
intermediate unless the energy dissipation is taken into 
account. The growth of a metastable domain is not spon- 
taneous but forced by the external field. The final state 
is determined merely by how many photons are absorbed 
to increase the ionicity. This result is consistent with the 
experimental finding in the neutral-to-ionic transition in- 
duced by intrachain charge-transfer photoexcitations. 21 ) 
By intramolecular excitations, on the other hand, a tran- 
sition from the neutral phase to the ionic phase with 
three-dimensionally ordered ferroelectric polarizations is 
shown to be induced, as clearly demonstrated by X-ray 
diffraction. 25 ) The reason for this difference is unclear 
yet. 



The qualitative difference described above between the 
two transitions is not due to the different inversion sym- 
metry between the dimerized ionic and the regular neu- 
tral states because the difference survives even if the 
electron-lattice coupling is turned off so that both ionic 
and neutral states have even parity with respect to the 
space inversion. The most plausible reason for the differ- 
ence is the fact that the ionic state is a Mott insulator 
caused by the electron correlation, whereas the neutral 
state is a band insulator caused by the band structure. 
In a Mott insulator, all the electrons are so correlated 
that any one of them cannot easily make the first move. 
However, above the threshold excitation density, they 
cannot tolerate the increased total energy. Once some 
electrons move, they trigger the collective motion. One 
electron is transferred after another like dominoes. The 
collective dynamics in the ionic phase are due to the 
electron-electron interaction that realizes the Mott insu- 
lator phase. They cannot be described by stochastic ap- 
proaches of refs. 4 and 5. In a band insulator, electrons 
move individually. The supplied energy is merely con- 
sumed to transfer electrons from the donor to acceptor 
molecules almost independently They do not strongly 
influence the motion of other electrons. 

Thus far, we have focused on the dynamics in purely 
one-dimensional systems. Although the short-time be- 
havior may not suffer from interchain interactions, they 
are important for the condition of coherence, which can 
be manifested by a double pulse. 26 ) The coherent mo- 
tion of the macroscopic neutral-ionic domain bound- 
ary 18 ;) - g possible only when interchain interactions are 
strong enough. Although the effect of interchain elas- 
tic couplings is considered in the context of confine- 
ment of a metastable domain, 17 ) interchain electron- 
electron interactions are much stronger. The intramolec- 
ular charge distribution in the TTF and CA molecules is 
calculated by the ab initio quantum chemical method. 27 ) 
The interchain electrostatic energies between neighbor- 
ing molecules are smaller than but comparable to the 
intrachain ones. Because the molecules arc tilted, the 
relative repulsion strengths are very different from the 
naive expectations based on the intermolecular distance 
measured in terms of the center of mass of each molecule. 
In the ionic phase, the interchain attractive interaction 
between the neighboring donor and acceptor molecules 
along the b' axis is larger than any interchain repulsive 
one between the donors or between the acceptors. The 
consequent electrostriction is theoretically shown 28 ) to be 
responsible for the pressure-temperature phase diagram 
of the TTF-CA complex containing the paraelectric ionic 
phase 29 ) and for the discontinuous contraction along the 
b axis. 30 ) 

Then, we add to the model (10) with (14) interchain 
electron-electron interactions , 3 1 ) 

#k!ter = X! (U±Sni d Sni^ + i + V±iSni d+1 5ni +1J ) , 

hi 

(16) 

where j is the chain index, and the donor-acceptor cou- 
pling strength Vj_i is slightly larger than the donor-donor 
or acceptor-acceptor one £/j_. 27 ) When the ionic phase 
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is photoexcited, the transition dynamics depend on the 
strengths of the interchain couplings. With weak inter- 
chain couplings, for instance, halves of the values for 
TTF-CA, the interchain correlation is very weak during 
the transition. A first neutral domain is easily created by 
a low density of photons. If the electric field is switched 
off immediately after this absorption, the residual chains 
remain ionic. By continuing the application of the elec- 
tric field, next domains are created in the neighboring 
chains. Many photons are finally needed to complete the 
transition into the neutral phase. With strong interchain 
couplings comparable to those in TTF-CA, the inter- 
chain correlation is strong during the transition. To cre- 
ate a first neutral domain, more photons need to be ab- 
sorbed. Once a domain is nucleated, nearby domains are 
almost simultaneously created in the neighboring chains 
to grow cooperatively. Their growth cannot be stopped 
even if the electric field is switched off immediately after 
the appearance of the first domain. As a consequence, 
the density of absorbed photons needed to complete the 
transition is lower than in the weak-coupling case. Be- 
cause of the strong interchain correlation, the growth of 
metastable domains coherently proceeds, which is consis- 
tent with the experimentally observed, coherent motion 
of the macroscopic neutral-ionic domain boundary. 18 ) 

3.5 Photoinduced phases in MX and MMX chains 

Besides the TTF-CA complex described above, many 
other molecular materials have been or are currently 
studied to show their own transition dynamics. Some 
of them show very asymmetric transitions between two 
different electronic phases. In a halogen-bridged nickel- 
chain compound, the ultrafast photoinduced Mott tran- 
sition from a charge-transfer insulator to a metal is 
observed. 32 ' Photoinduced optical responses of one- 
dimensional half-filled Hubbard models with and with- 
out alternating potentials are studied by the exact diag- 
onalization method to show a large shift of the spectral 
weight to low energies including the Drude component 
in the Mott insulator phase, which are similar to the re- 
sponse caused by chemical doping. 33 ' 

In halogen-bridged binuclear platinum-chain com- 
pounds, R4[Pt2(pop)4l]nH 2 (pop=P 2 05ll2~) with 
cation R, the transition from a charge density wave 
(CDW) to a charge-polarization (CP) state is photoin- 
duced much more easily than its inverse transition. 34 ' 
The electronic phases of these compounds are controlled 
by the distance between the neighboring binuclear units 
diviXM both when the counterion R is varied and when 
the applied pressure is varied. The nonmagnetic CDW 
phase is observed for small (JmxMj while the paramag- 
netic CP phase for large ^mxm ■ The charge density is dis- 
proportioned among binuclear units in the CDW phase, 
and within each binuclear unit in the CP phase (Fig. 8) . 
In both phases, the iodine ion approaches the platinum 
ion with less electron density. The dependences of the 
electronic phase and its physical properties on <iMXM 
are well explained by perturbation theories from the 
strong-coupling limit as well as the exact diagonalization 
method. 35 ' The mechanisms of the photoinduced phase 
transitions are clarified by taking the charge transfer 
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Fig. 8. Schematic electronic and lattice structures in MMX 
chains. 

processes explicitly into account 36 ' in a one-dimensional 
two-band three-quarter-filled Peierls-Hubbard model, 

H = - t M M X^ C 2Z-1 i( t C 2Z,<t + n ' C ') 

1,(7 

- iffiM^fc^/a+Lu + h.c.) 

+ f3}^qi(n 2 i+i - n%) + UjA/,ni^m,i 
i i 

+ K MX J2qi+(M/2)J2<8 > ( 17 ) 

I I 
where c, creates an electron with spin a at site I, 

ni.a = c i a c l,a, and n; = n i.a- The binuclear unit con- 
tains two M (M=Pt) sites, 21—1 and 21. The displacement 
of the X (X=I) ion between the two M sites 21 and 21+1, 
relative to that in the undistorted structure, is denoted 
by qi . The transfer integral within the unit is denoted by 
*mm, while that between the neighboring units through 
the X p z orbital by £mxm- The energy level of the M d z 2 
orbital depends on q with linear coefficient (3. The on-site 
repulsion strength is denoted by Uyi, the elastic constant 
for the MX bond by -KjyiX) and the mass of the X ion by 
M. With increasing pressure corresponding to decreasing 
c?mxm ■ P increases because it originates from the electro- 
static energy between the positive M and the negative X 
ions. This is why the discontinuous transition from the 
CP to the CDW phases is induced by pressure. 34 ' 

It is true that the dependence of the threshold excita- 
tion density in the photoinduced CDW-to-CP transition 
on the pressure, which controls the relative stability of 
these phases, can be explained qualitatively with their di- 
abatic potentials, i.e., in a stochastic manner. However, 
the asymmetry in these transitions requires explanation 
based on the explicit charge transfer processes (Fig. 9): in 
the CP state, low-energy charge-transfer processes occur 
only within a binuclear unit, and consequently they do 
not lead to the CDW state. Even if the CP phase is irra- 
diated by higher-energy photons, charge transfer among 
the binuclear units takes place incoherently. Since the in- 
tcrunit charge transfer is energetically unfavorable, the 
CDW domains never proliferate. Furthermore, the tran- 
sition amplitude connecting the degenerate CDW phases 
with opposite polarizations is much smaller than those 
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Fig. 9. (Color online) Schematic ground and low-energy photoex- 
cited states in MMX chains. 36 ' In the CDW phase, the charge is 
transferred to the neighboring unit to form a transient CP do- 
main. In the CP phase, the charge is transferred within the unit 
to merely reverse the polarization. 



connecting the degenerate CP phases. The coherence re- 
covery is achieved with much more difficulty in the CDW 
phase. 

4. Concluding Remarks 

We have summarized the history of theories for pho- 
toinduced phase transitions, focusing on the development 
from stochastic to coherent transition dynamics. Both 
dynamics are characterized by cooperativity. Coherent 
dynamics are direct consequences of interactions among 
electrons and/or between electrons and lattice displace- 
ments. Photoinduced transition dynamics between differ- 
ent electronic phases of different materials will be stud- 
ied further in appropriate models for understanding each 
mechanism and for exploring new possibilities. Current 
and future issues to be solved will include differences be- 
tween photoinduced and thermally induced phases, ori- 
gins of different nonlinear characters, coherence originat- 
ing from phonons and/or that from electronic excita- 
tions, relative dynamics among spins, charge and lattice, 
relations to quantum phase transitions or quantum na- 
tures of phonons, and how to control the dynamics by 
tuning the laser pulse. 
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